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ABSTRACT 


Low density nozzles or large expansion ratio nozzles used in space experience 
rarefaction effects near their exit in the form of velocity slip and temperature jump 
at the walls. In addition, the boundary layers become very thick and there is 
a very strong viscous/ inviscid interaction. For these reasons no existing design 
technique has been found to accurately predict the nozzle flow properties up to the 
nozzle exit. The objective of this investigation was to examine the slip boundary 
conditions and formulate them in a form appropriate for use with a full Navier 
Stokes numerical code. The viscous/ inviscid interaction would automatically be 
accounted for by using a compressible Navier Stokes code. Through examination of 
the interaction of molecules with solid surfaces a model for the distribution function 
of the reflected molecules has been determined and this distribution function has 
been used to develop a new slip boundary condition that can be shown to yield 
more realistic surface boundary conditions. 
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INTRODUCTION 


As the density in a flow field diminishes the continuum description of that 
flow field begins to break down, first at the wall, and later in the interior of the 
flow field. The relevant parameter which describes the breakdown of the continuum 
description is the Knudsen number, defined as the ratio of a characteristic molecular 
mean free path to a relevant length that is descriptive of the flow field. The Knudsen 
number can be related to the Mach number and the Reynolds number, the simplest 
argument indicating that 



The relevant mean free path and length scale depend upon the physical problem 
(see ref. 1-6). 

If the Knudsen number is not considerably less than one, the flow is no longer a 
continuum. The most significant feature of the flow when there is slight rarefaction 
is a change in the boundary conditions at the wall. This regime is characterized 
by the fact that the rarefaction effects (i.e., changed wall boundary conditions) are 
small corrections to the continuum flow behavior. However, viscous effects become 
more important as Kn increases. Then the skin friction becomes comparable to the 
wall pressure and the heat transfer coefficient increases. The displacement effect of 
the thick boundary layer changes the surface pressure from its inviscid value. These 
are the dominant fluid mechanical effects in the slip flow regime, rather than the 
relatively small corrections introduced by the modified wall boundary conditions. 


The transition regime occurs when the Knudsen number is of the order of 
one. In this regime molecular collisions with the surface and with the free-stream 
molecules are of more or less equal importance. The flow is no longer a slight 
correction to the continuum flow. At still larger values of the Knudsen number 
only collisions between the incoming molecules and the body surface need to be 
considered. This latter regime is called free-molecule. The delineation between the 
various regimes, i.e., continuum, slip, transition and free-molecule, is not sharp and 
depends upon the Mach number and the body geometry. References 1-6 should be 
consulted for more details. 


Rarefied gas dynamics is the field that describes these flow regimes. NASA 
obviously has many applications for rarefied gas dynamics. These include the de- 
sign of low density nozzles or large expansion ratio nozzles for use in space, the 
determination of contamination to spacecraft solar panels, optical sensor windows 
or temperature control surfaces, the estimation of drag and torque on satellites or 
the space station and of the heating and aerodynamics of vehicles entering or skip- 
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ping along the edge of the atmosphere. The present investigation concentrated on 
those applications that required slip flow boundary conditions. 

The objectives of the present work are presented in the following section. This 
is followed by a brief review of kinetic theory, which is the tool used to solve rarefied 
gas dynamics problems, and of gas-surface interactions. The latter is required to un- 
derstand the new formulation of the slip boundary condition that was accomplished 
for this study. The new slip condition is then presented and compared with previous 
slip boundary conditions. Evidence for the superiority of the present formulation is 
provided, but only incompletely. Finally, recommendations for additional work are 
made. 
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OBJECTIVES 


NASA has a continuing need to improve nozzle contour design. Many of 
these nozzles are used under low density conditions but no design technique that is 
presently available can predict the flow field near the exit plane under these condi- 
tions. Accurate exit plane prediction is required to estimate thrust, as an input to 
contamination prediction codes and to accurately determine the back pressure on 
end plates of ultra-low thrust satellite control rockets. 

With these applications in mind, the objective of this investigation was to 
devise an accurate means for the design of low density nozzles. To accommodate the 
strong viscous/inviscid interaction that will exist when the boundary layer becomes 
very thick near the nozzle exit, it will be necessary to use a full compressible Navier 
Stokes code and to add slip boundary conditions when the Knudsen number becomes 
sufficiently large. A second objective was to utilize recent gas surface interaction 
studies to develop a model for the distribution function of the molecules reflected 
from the surface that could be used to determine a slip boundary condition that 
can more accurately predict phenomena that are observed in slip flow, including 
some anomalous effects which have been observed in capillary flow. 
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PRINCIPLES OF KINETIC THEORY 


Kinetic theory is the discipline that uses the molecular description of gases to 
determine the flow properties of a gas that is not in equilibrium. Because the number 
of molecules in any finite but small volume is very large, a statistical approach can 
be taken. The fundamental entity in kinetic theory is the distribution function, /, 
which is a function of position (x,y,z) , velocity ( u,v,w ) and time, t, such that 

f(x,y,z,u,v,w,t)dxdydzdudvdw (l) 


is the number of particles in the volume dxdy dz about position (x, y, z) with velocity 
in the range dudvdw about ( u,v,w ) at time t. If the distribution function is known 
then all flow properties can be obtained as integrals (moments) of the distribution 
function. For example, letting dr = dxdydz and du = dudvdui 


n(r,t) 


OO 

/// /(r ’“' 


t)du, number density 


( 2 ) 


OO 

u = — J j J f (r.u,t)du, mean (gas dynamics) velocity 


( 3 ) 


OO 

"--/// U t U 3 f{r,u,t)du 


( 4 ) 


P l} is the flux of i-th momentum in the j-th direction. U is the peculiar or random 
velocity, defined as 

U = u -u . (5) 

Other properties can be similarly defined; see reference 8. 

To proceed further a molecular model must be developed. The simplest model 
is that of molecules composed of rigid elastic spheres of diameter a which interact 
with one another only upon impact This is the model that will be assumed in this 
work while more realistic models are described in reference 8. Classical mechanics 
can then be used to show that the distribution function varies with time according 
to the equation 


£1 dj_ 
Dt dt 


+ u • 




-C<~> = J(/,/ t ) 


(6) 


IX -6 



where C dr du dt is the total number of molecules having velocities in the range 
u to u + du which are lost in dr du in time dt due to collisions with molecules 
having other velocities and C dr du dt is the total number of molecules gained 
to the range a to a + du in dr du in time dt due to collisions between molecules 
having those initial velocities that lead to a final velocity of one of the molecules in 
the range u tou + du. The specific form of the quantities C and is given 

in reference 8. It has been assumed that the molecules are not acted upon by an 
external force field. 

The above equation is called the Boltzmann equation and is considered to be 
the fundamental equation of rarefied gas dynamics. Many assumptions must be 
made to derive this equation, the two most important being that the equation as 
written applies only to a pure gas with no internal degrees of freedom (monatomic 
molecule). The equation can be generalized to eliminate these limitations (ref. 9) 
but even in this simplest case illustrated the functions and are so complex 
that the equation can be solved only for a few special circumstances. 

Moments of the Boltzmann equation, i.e., equations of the form 

ill - -/(/.A)} = °. (7) 

— OO 


play a special role in kinetic theory. These are the Maxwell equations of change. 
Since mass, momentum and kinetic energy are conserved in each molecular collision 
(monatomic molecules) the moments of the collision integral for n = 0, 1 and 2 are 
zero. The moments of the left hand convective term lead to the familiar continuity, 
momentum and energy equations of continuum mechanics. Each equation, however, 
contains a higher order moment and the system is not closed. The system can be 
closed for systems close to equilbrium in the limit as Kn goes to zero, leading to 
the Chapman-Enskog expressions for the stress tensor and heat flux vector. 

A useful exact solution of the Boltzmann equation exists for a gas in a state 
of local equilibrium, called a locally Maxwellian distribution (ref. 8). It can be 
written as 



where n and u depend upon r and t and T depends upon t. The collision integral 
is zero for the locally Maxwellian distribution, i.e., •/(/o/o 1 ) = 0, and the stress 
tensor and heat flux vector are zero. The conservation equations become the Euler 
equations of fluid mechanics. 
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A useful property of a gas is the mean free path, which is defined as the mean 
distance that a molecule travels between collisions. This concept has meaning only 
for a molecular model of hard elastic spheres. If the gas is described by a Maxwellian 
distribution then the mean free path is given by the expression 

A = — - . 

\/2 n 7 T o~ 

More complicated situations are discussed in reference 8. It was assumed in the 
Introduction that a single mean free path described the flow but that is usually not 
the case. 

The purpose of this brief introduction has been to indicate that boundary con- 
ditions must be provided on the distribution function itself to solve problems in 
rarefied gas dynamics. For problems of interest the most difficult boundary condi- 
tion to specify is that at the wall. The interaction of molecules with solid surfaces 
will be reviewed in the next section and it will be shown that the experimental re- 
sults can be described by a particularly useful model of the distribution function of 
the reemitted molecules. This model will then be used to develop a new boundary 
condition in the following section and this will yield a new slip regime boundary 
condition. 
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GAS-SURFACE INTERACTIONS 


The study of gas-surface interactions provides the boundary conditions for the 
calculation of low density flows past wall boundaries. Our understanding of the 
molecule/ wall interaction has evolved with time and is still presently incomplete. 
The range of relative energies between the molecule and the surface for applications 
of interest is very broad. Satellites operate in a free-molecule environment with 
a relative velocity between the satellite and the surrounding gas of around 26,000 
ft/s (approximately 8000 m/s). This translates into a relative energy of 9.leV for 
a nitrogen molecule or 5.2eV for an oxygen atom. The flow in low thrust or low 
density nozzles offer examples of the other extremes. The flow near the exit of a 
low density nozzle which is operating in space is in the slip regime on the wall. 
In this case the energy of the gas molecule relative to the wall is considerably 
lower, frequently below O.leU, even though the wall and gas may be at a high 
temperature (the average molecular energy at 3000°K is 0.39eV) This wide range 
of relative energy of interest has caused experimental problems with the study of 
the interaction of gases and surfaces because so many physical processes can occur 
over such a broad energy range (ref. 10) and because of the difficulty of producing 
relative energies greater than 0.3eU; in fact, the satellite energy range has been 
reached only comparatively recently in laboratory experiments. 

The study of gas-surface interactions has a long history that can only be high- 
lighted here. The purpose of the discussion will be to justify the form of the distri- 
bution function of the reemitted molecules that will be used to obtain the slip flow 
boundary conditions 

Early experiments indicated that there was some accommodation between the 
incident molecules and the wall properties. If molecules reflect specularly then there 
is no accommodation to the wall conditions and Maxwell reasoned that a portion 
must be reemitted diffusely after having come to complete accommodation with 
the wall. He introduced the reflection coefficient, /, which is the fraction of the 
molecules which reflect diffusely while the fraction (I — /) reflect specularly (ref. 
II), Later it was recognized that different properties of the molecule-wall collision 
accommodate to a different degree and accommodation coefficients were introduced. 
The accommodation coefficients for normal momentum, tangential momentum and 
energy are defined by the equations (refs. 1 and 3). 


a 

a 


Pi -P 


Pi ~ Pt 


— = normal accommodation coefficient 


r t - r 


- = tangential accommodation coefficient 


e x — e 


— = thermal accommodation coefficient 


( 10 ) 

( 11 ) 

( 12 ) 


IX-9 



where p t , r t and e t are incident fluxes of normal and tangential momentum and 
energy, p r ,r r and e r are the reemitted fluxes of these quantities from the surface 
while p w and e w ( t w = 0) are the fluxes that would be reemitted by a gas in complete 
Maxwellian equilibrium with the surface. These expressions have significance only 
for conditions close to equilibrium. 

Almost all rarefied gas dynamic problems have been solved using these accom- 
modation coefficients with values that have been obtained from laboratory mea- 
surements. The measurements that correspond closest to the defined quantity are 
those of the thermal accommodation coefficient (refs. 12 to 14). Recently o and o ' 
have been measured by measuring the forces exerted on a small plane surface, that 
is connected to a vacuum microbalance, by a stream of molecules (molecular beam) 
of known energy. Both of these quantities have been measured as a function of the 
energy of the incident molecules and of the angle between the surface normal and 
the incident molecular beam direction (E l ,0 l ). 

The low energy and satellite energy behavior of o and o' have been found 
to be quite different (refs. 15-18). The low energy behavior was examined very 
completely by Seidl and Steinheil (ref. 17). In the "as received "condition plus 
an ambient temperature desorption, the surface of all materials was found to be 
covered with adsorbed gases, primarily water vapor and hydrocarbons However, o 
was found to depend greatly upon the bulk material even though the surface layers 
were the same. The tangential accommodation coefficient was greatly diminished 
for all values of 9 t by removing the adsorbed gases if the surface was atomically 
smooth but o was greater than one at small 0 t for atomically rough surfaces due to 
back scattering. 

The tangential accommodation coefficient was found to fall with increasing 0 t 
at low energies (ref. 17) but the opposite behavior was observed at satellite energies 
(ref. 18). Also, the integrated values of o' (over all 6 t ) decreased from their low 
energy value near unity to a minimum at an energy somewhere below satellite energy 
and then increased to values comparable to the low energy value at higher energies 
(refs. 10, 15 and 16). 

The measurements of a and o' at satellite energies by Knechtel and Pitts (ref. 
18) are of interest. They measured both coefficients as a function of 0 t and E t for 
nitrogen on aluminum and then applied them to determine the torque and drag 
on a surface element in free molecule flow at satellite energy (E t = 9eV'). A very 
large effect was shown to occur for the torque when compared to the assumption 
of full accommodation (see Fig 16 in ref. 15) The drag was also calculated for a 
spherical satellite as a function of E t . For free-molecule flow in the limit 5 — » oo 
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it can be shown that the drag coefficient on a sphere reduces to (ref. 1) 

C D — > (2 - a + a). (13) 

In either the limit of full accommodation [a = o' = l) or no accommodation 
(cr — <j f = 0) this gives the value Co ~ 2, and in general Co > 2. However, 
integration of the forces using the measured values of a and a 9 yielded values of 
Cq less that 2 at E t less than 40^. This significant finding will be compared with 
calculations that do not use accommodation coefficients later. 

The momentum accommodation coefficients were historically defined only for 
free molecule flow. But they are also important for the understanding of gas- 
surface boundary conditions in transition and slip flow (ref. 19). The velocity slip 
and the onset of rarefaction effects can be significantly influenced by the degree of 
accommodation (ref. 20). 

The accommodation coefficients only allow integrated properties of the reflected 
molecules to be determined in terms of the incident properties (fluxes are integrals 
over the distribution function - see the next section). However, what is desired is a 
joint distribution function (ref. 1) 


(14) 

which is the probability of a molecule incident upon a wall with velocity u t being 
reflected with velocity u r . This distribution function would be a function of nu- 
merous parameters, such as the type of surface material, the microscopic surface 
state, the wall temperature and the presence or absence of adsorbed gases. Present 
knowledge of the gas-surface interaction process is insufficient to yield knowledge 
of the joint distribution function. However, molecular beam scattering experiments 
have given considerable information about the interaction process. The most useful 
experiments are those that impinge a surface of known state with a monoenergetic 
molecular beam at various 0 t and measure the reflected number density and velocity 
distribution as a function of reflected angle (ref. 21). These experiments have not 
yet produced enough information to give g, however they have assisted in the devel- 
opment of a parametric representation of the distribution function of the reflected 
molecules. This representation will be used in the next section to develop a new 
slip condition while the experiments themselves will be examined to estimate the 
values of the parameters in the representation. 

Accommodation coefficients have significance only for conditions close to equi- 
librium (ref. 13) Frequently flows of interest are far from equilibrium and the 
results supplied by the accommodation coefficients are not adequate (ref. 21). The 
molecular beam simulates conditions far from equilibrium and should thus model 
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such conditions better than standard accommodation coefficients. 

Considerable progress has been made toward determining the joint probability 
distribution. Recent molecular beam scattering experiments are described in refer- 
ences 22 to 29 while the earlier experiments are summarized in reference 13. Only 
certain aspects of these measurements are relevant to the present discussion and 
they can be summarized as follows. 

1. Maxwellian diffusive reflection, with the molecules accommodated to the 
wall temperature, occurs at low beam energies if the energy of adsorption of the 
molecules is high or if the surface is covered with adsorbed gases. 

2. Lobal rather than diffusive (cosine) reflection with incomplete thermal 
accommodation is observed at low energies if the surfaces are properly degassed 
or if the incident beam energy is increased such that the incident kinetic energy is 
greater than the wall temperature. 

3. Lobal reflected distributions are always observed at higher energies whether 
the surfaces are degassed or not and from polyecrystalline or single-crystal materials. 

4. The peak in the lobal distribution is at an angle 6 m from the surface normal 
which is less than the corresponding incident angle of the molecular beam, 0,, at low 
energies. As the energy of the incident beam increases (for fixed 0 t ), Q m increases 
so that the lobe peak approaches the surface tangent, but at some energy less than 
lOeV, d m reaches a maximum value and decreases toward the specular angle with 
further increases in the beam energy. This pattern is approximately the same for 
all e t . 

5. The speed distribution of the reflected molecules is analyzed by fitting 
the data to a " drifting "Maxwellian with parameters U r and T r . U r is the mean 
velocity of the distribution. U r ,U t and the surface normal are in the same plane. 
T r is a measure of the spread of the distribution and is not usually equal to the 
wall temperature T w . The parameters are usually combined to give a speed ratio 
S T = U r /(2kT r /m) 1 / 2 . Earlier experiments (ref. 30) indicated that a one-parameter 
distribution was insufficient to characterize the reflected speed distribution but most 
data can be adequately characterized using these two parameters. 

6. At low energies, as 0 r — * 90° (tangential reflection), S r increases and T r 
decreases for given S t and 0, (refs. 26 and 28). 

Some years ago Nocilla (ref. 31) suggested that the reflected molecular distri- 
bution could be assumed to have the form of a " drifting "Maxwellian of the 


form 


/»» 


n r 


/ m 
\ 2nkT r 



-m(u-U r ? ] 
2kT r j 


(15) 


where u is the molecular velocity, U r is the macroscopic drift velocity attributed 
to the outgoing stream which is assumed to lie in the plane in the incident velocity 
and the surface normal, T r is the effective temperature of the distribution (not nec- 
essarily the wall temperature T w ) and n r must be found by balancing the incoming 
and outgoing number fluxes. The previous discussion indicates that this form of 
distribution function has been found to describe most of the molecular beam data. 


Nocilla made a comparison between some of the scant angular flux data that 
were available at the time and determined a best fit using this distribution. The 
fit produced the parameters 6 r and S r (the distributions were normalized to a unit 
value of (2 kT r /m) 1 ^ 2 ). The only measurements that existed were for the scattering 
of low energy nitrogen from lithium fluoride crystals. Polished crystal surfaces tests 
produced 6 r > 90° and small values of S r . Cleaved crystal surface tests produced 
0 r slightly greater than 0, and S r of the order of 1 0; an exception occurred for 
0 t approaching 90°, for which 0 T > 90°. All the measurements were adequately 
described by this function except the tests with 0, near the surface tangent. In this 
case the data was adequately represented except in the quadrant that contained 
U t , where the distribution underpredicted the scattering flux. In this case a greater 
number of molecules were rejected toward the source than the distribution function 
predicted. 

This distribution function was used by Hurlbut and Sherman (ref. 32) to 
determine the drag on a spherical satellite. The value of T r was determined from an 
angle-dependent thermal accommodation coefficient. Both S r and a were assumed 
to vary linearly with 0 t ; U T was assumed to be directed at 9 r , which was taken 
to be equal to 0 t . The results yielded the classical diffusive or specular reflection 
results for appropriate values of the parameters but notable among the results was 
a prediction that Cq could be less that 2 0 for S T — + oo for some values of the 
parameters. Recall that this result was also predicted by Knechtel and Pitts (ref. 
18). 

The Nocilla distributions will be assumed to be an accurate representation of 
the reflected distribution function for slip flow and a new slip boundary condition 
will be presented. The problem of the estimation of the parameters U r ,S r and 0 r , 
will still remain, but the experimental results described in this section can be used 
as a guide. 
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SLIP WALL BOUNDARY CONDITIONS 


Slip flow occurs when there is only slight rarefaction and the departure from the 
continuum behavior described by the Chapman-Enskog solution is small. The most 
rigorous way of solving this problem near the wall would be to solve the Boltzmann 
equation using an appropriate kinetic boundary condition (ref. 8). Instead an 
approximation due to Patterson will be used (refs. 6 and 33). 

The molecular motion near a wall results from the interaction of the incident 
and reflected molecules. The incident molecules within a mean free path of the 
wall have a motion consistent with the external flow while the reflected molecules 
have come to some partial accommodation to the state of the wall, as discussed 
previously. The interaction of the incident and reflected molecules results in a 
deviation in the distribution of molecular velocities from a local Maxwellian. This 
distribution function next to the wall must be obtained by solving the Boltzmann 
equation, which is a very difficult task. 

A few mean free paths from the surface the gradients of the macroscopic velocity 
and temperature can be considered to be constant. These gradients cease to be 
constant only in the immediate vicinity of the wall where the collisions between the 
incident and reflected molecules have their effect (see figure 1). The extrapolation of 
the gradient to the surface yields the slip velocity or extrapolated wall temperature. 
Note that this does not yield either the gas velocity or the gas temperature at the 
surface. But as far as the velocity and temperature fields in the bulk of the flow 
are concerned, these quantities can be used as boundary conditions for the Navier 
Stokes equations (second term of the Chapman-Enskog expansion) to yield accurate 
solutions everywhere except within a few mean free path lengths of the wall. 

Patterson assumed that the gas near the wall was not in equilibrium but was 
close to equilibrium. The distribution function was assumed to be of the form 


/. = /o[i + *Wi 


H i 


/ m \i/2 

t m \ l /2 


(*f) u ' tHi = 

51 

£ 

ii 

(*t) 


1/2 


U z 


fo is the local Maxwellian. 


(16) 


— - ^ H 2 

f 0 dU x dUy dU z = (2tt)- 3/2 e ' * ‘ dHidH 2 dH-i (17) 

3 3 3 3 3 3 

F(ff) = i ff, + g £ £ < l8 > 

1 11 111 
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where the order of the subscripts is immaterial, e.g., a 12 = a 2 i and a 122 = a 212 = 
a 22 i. n, T and U 2 axe defined by integrals of /o. 

Certain conditions must be satisfied by the coefficients, namely (ref. 33). 

°11 + 0,22 + a 33 = 0 (19) 

and 

2a, + o, 1 1 + a t 22 + o, 33 = 0 (f = 1, 2 or 3) (20) 

Patterson demonstrated how to relate the remaining coefficients to integrals involv- 
ing properties of the intermolecular collisions which were evaluated for hard elastic 
spheres. 

Now consider the interaction of the incoming and reflected molecular streams 
in the vicinity of the wall. The motion of the gas on the surface depends on the 
transfer of mass, momentum and energy by the gas molecules to and from the 
surface. Consider figure 2 with the axis fixed relative to the body. The volume S in 
the flowing gas contains gas molecules with mean velocity components u 3 , 0, uj a and 
a random motion described by /,. The reflecting molecules are assumed to act as 
though they emerged from an adjacent volume W with a Nocilla distribution with 
U r lying in the (x — y) plane. The molecular motions in the volumes S and W are 
mutually compatible if the total mass, momentum and energy transferred upward 
by the molecules from W to S are respectively equal to the total mass, momentum 
and energy which would be transported upward if the molecules in W moved with 
the same mass velocity and random motion as the molecules in S. 

Compatibility of the mass flux yields 


OO OO OO 


OO OO OO 


Iff v f t (T a )du dv dw = mn r /// v f r (T r )du dv dw 


— oo 0 — oo 


-oo 0 — oo 


( kT r \ [ m(u-U r ) 2 ' 

= “ P [ 2*7V — J 

Integration and the use of the conditions on the a coefficients yields 

mn a , mn r , . 


x(o r ) = exp (—of) + \/no r [l + erf(o r )j 


2 kT r 


( 22 ) 

(23) 

(24) 

(25) 
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This result can be written as 


n r [Ta |1+ a22/2] 3 , , 

* V TV xM <26) 

In the hard sphere approximation used by Patterson, 

5 [ A f du dw dv\] 

+ *-**)]. < 27 > 

The balance of the transport of the tangential component of momentum at the 
interface yields an equation for the slip velocity. 


OO OO OO 


OO OO OO 


' TI 3 Iff uvf l (T 3 )du dv dw = m n w ffl uv f r (T r )dudv dw (28) 


-00 0 — 00 


(2tt)3/2(2^) 


-00 0 —00 


7T \fl/K dy 2 + &122 f + 

J J 3 


a 2 3 /2 v /2^d 
n r m U? . 


t(2 + a 2 2)l a (29) 


x(o r ) sin 6 r 


~ 2yfn S T v ' 

Using the previous result for n r fn s and dropping higher order terms yields 


kT r 

u 3 = —2 S r \ — - sin 6 r — 


V 27 r 012 4- 


5 A / du dv 
12 4 [f \ dy dx 


15 / tt\ 1/2 A dT 
0,22 ~ 16 V 2/ T dx 


W 3 = - 


032 + O 322 


5 A / diu dv 


0-32 4 6 \ dy d2 / 

15 / 7r \ l /2 A dT 
1322 - _ 16 \2/ T dl 
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The balance of the transport of the normal component of momentum at the interface 
yields a relation between p 3 and p w . 


oo oo oo 


oo oo oo 


m 


-III v 2 fi(T 3 )du dv dw = mn r Iff v 2 f r (T r )dudv dw (36) 


— oo 0 — oo 

Integration yields 


— OO 0 — OO 


Ei 

2 


1 + O22 + 


3 \Z 2 tF . 


= ^ { a '*M 


+ —■ (1 + erf cr r ) | 


(37) 


where the perfect gas law 


p = n kT 

has been used. This can be written as 


El 

Pr 

0 2 22 


1 — 022 — 


2o 2 


22 


3 \/ 27 r J 

_45 /ttN 1 / 2 A dT 
~~16\2/ Tdy 


1 + erf Or + — 7= °r X (<?r) 
V 71- 


( 38 ) 

( 39 ) 


The temperature jump between the gas and the wall can be determined by the flow 
of translational energy of the molecules across the interface. 


OO OO OO 


i m n 3 Iff v(u 2 + v 2 + i v 2 )f t (T 3 ) du dv dw 


— oo 0 — oo 


oo oo oo 


(40) 


= ^mn r /// v(u 2 4- v 2 + w 2 ) f r (T r )du dv dw 


— oo 0 — oo 

Integration yields the following result 
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— mn r 


(S r 2 + 2 )x{&r) H — — cr r (l -j- erf 


' m \2nm J ^ r "" ' 2 ' 

Only keeping higher order terms, assuming S r is small, gives 


T a f 1 1 pw 1 f / S?\ [tt o r (1 + erf o r ) ) . , 

T r = [ t -4* + dl«l{{ 1+ fHi XM 1 <42) 

75 /7 rs 1 / 2 A ar 

02 ~ 32 \ 2/ ^ 

These results can be compared with the results of other investigators by writing the 
expressions in the form used by Kogan (ref. 34.). He wrote, for diffuse reflection 

/x / 2 du 2pi 1 dT 

U * = A n V mkT dy + nm T dx 

K / 7rm dT 

AT* = C */ — (45) 

2A:n V 2A;T dy V ' 

where /x is the coefficient of viscosity and K is the coefficient of thermal conductivity. 
The magnitude of the coefficients A, B and C depend upon the equation used for 
their determination and the molecular model assumed. 

The present results can be written in this form by using the transport coefficient 
relations derived by Patterson for a hard sphere molecular model, namely, 

x 1 r.™ / L-'T \ 1 / 2 


•if! ™ 


A = 

1 


5m 


yf2nno 2 

/x = 

16 <7 2 

_ 

/ 8 /tT \ 1/2 


15 

(7 = 

\7r m / 

K = 

■J A* ' 
4 i 


Then, the results can be compared as follows: 

1. Patterson, hard spheres, diffuse reflection. 


„ M / 2 ( du dv 

0.8862 - \l — + — 

n V mkT \ dy dx 


A 

m 


„ 2/x 1 dT 

0.3750 — 

nm T dx 


K j 7 r m dT /x / dtx du; dv 

1 + ^ - 0 1667 pU + 57 


or A = 0.8862, fi = 0.3750 and C = 1.000. 
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2. Present results, using Nocilla reflection model and hard spheres. 


kT r I u 

u 9 = —2 S r \j sin 0 r + { 0.8862 
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2 / du dv 

I h — 


n V mkT \dy dx 


Ts 

T w 


=1 


1 + 


K 


+0.3750 


2 n 


nm 


nm dT y 

0.1667 ~ 


2 km V 2kT dy 


i ar) 



T dx J 

> 
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(du 

dw 

_ o 


dz 

dy 

R °T 

(1 + 

erfcr r ) 


X(^r) 


( 49 ) 


(50) 


3. The traditional approach using accommodation coefficients and Patterson’s 
technique is given by Shidlovskiy (ref. 6). 0 = percentage of diffusely reflected 
molecules and a is the thermal accommodation coefficient. 


u s 


0.8862 


2 — 0\ a j 2 f du dv\ 2y 1 dT ( 

n\'^f U + ^) +0 ' 3750 ™ra^ (51) 


I± 

T w 



2 — 0 : K 
a. 2 kn 


nm dT 
2kT dy 


- 0 1667 - 

P 


( du dw dv \ 1 
\ dx dz ^ dy ) J 


(52) 


These results are identical to those of Chrusciel, et al. (ref. 18) except that their 
results were nondimensionalized. 


Kogan summarized other determinations of the coefficients A,B and C. This 
is repeated in Table I. 

The new slip formulation (equations (38), (39), (27), (49) and (50)) contains 
the parameters S r ,U r or T r and 0 r that must be specified before a problem can be 
solved. The parameters used will depend upon the problem being solved. That is, 
flow in a nozzle has the free stream flow tangent to the surface at all locations and 
it would be expected that 6 r would be close to 90° or greater, S r would be large and 
T r would be less than T w . On the other hand, for the flow over a blunt body, the 
parameters would be a function of 0 t and various forms could be assumed, much as 
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Hurlbut and Sherman did (ref. 32) for sphere drag, and ranges of the values could 
be examined for comparison with experiments. Notice that the form of the bound- 
ary conditions is somewhat different from those that include the accommodation 
coefficients. The slip velocity contains an subtractive term which is of the order of 
U r . The scattering experiments have shown that U r is never larger than the mean 
thermal velocity of the reemitted gas. For tangential flow over a surface this is 
small compared to the mean thermal velocity of the gas at the wall temperature. 
Since this term is negative it diminishes the magnitude of the usual slip boundary 
condition. This was also found by Chrusciel, et al. (refs. 21 and 35) using Monte 
Carlo methods to solve the flow over a blunt-nosed body in hypersonic flow. The 
temperature jump expression is multiplied by an expression that goes between the 
values 0.82 at o r = —0.20 to 1.44 for o r — 0.50. Values less than 1.0 would be 
required to match the results of Chrusciel, et al. 


Table I 


Coefficient 

Value 

Method of Determination 

A 

1.012 

Krook equation 

B 

0.383 

Krook equation 

C 

1.173 

Krook equation 

A 

1.047 

Boltzmann equation 
Maxwell molecules 

A 

1.015 


B 

0.383 


B 

0.402 


A 

1.103 

Boltzmann equation, 
hard spheres 

B 

0.4456 

\ \ 

B 

0.3292 

\ \ 

C 

1.113 


A 

0.8862 

Present results 

B 

0.3750 

Present results 

C 

1.000 

Present results 
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CONCLUSIONS AND RECOMMENDATIONS 


A new form of the slip flow boundary conditions has been developed which 
appears to have characteristics that can improve the agreement between the compu- 
tations and measurements. However, the boundary conditions have been obtained 
in an approximate manner and they contain unknown parameters. Therefore, the 
following steps should be taken to complete the task of obtaining useful slip flow 
boundary conditions. 

1. The analysis of the previous section should be repeated using the Chapman- 
Enskog expansion rather than the approximate solution of Patterson. 

2. Since the Chapman-Enskog expansion assumes that the flow is close to 
equilibrium, the analysis should be repeated by solving the Boltzmann equation 
close to the wall. It is not anticipated that this could be done exactly but a series 
expansion would be required. The boundary condition would then apply for flows 
arbitrarily far from equilibrium (ref. 21). 

3. The boundary condition should be used on a simple flow, such as Poiseville 
flow through a circular capillary. This problem is well understood theoretically (ref. 
36) and has an extensive data base (refs. 37 and 38). Therefore, it would be a good 
test case for the new boundary conditions and would help develop a rational for 
determining the parameters in the equations. 

4. Once the parameters have been obtained then the equations can be applied 
to the calculation of low density nozzle flows. For low throat Reynolds numbers 
(Re* < 100) the interaction of the wall viscous region with the nozzle core is so 
great that only the full Navier Stokes equations can be used to adequately solve the 
flow field. Therefore, an axisymmetric Navier Stokes code, such as that of Lombard 
(ref. 39), will be used with the new slip boundary conditions to calculate nozzle 
flow fields and compare the properties with any existing measurements that can be 
found (ref. 38). 
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